function Fig2_c_2()
Fcolor='k';
Bcolor='w';
Fsize=18;
Bsize=14;

%-1/2
data{1} = CS1_0_extraction_2D_Phi(46); %3T
data{2} = CS1_0_extraction_2D_Phi(50); %4T
data{3} = CS1_0_extraction_2D_Phi(44); %5T
data{4} = CS1_0_extraction_2D_Phi(48); %6T

%-3/2
% data{1} = data_extract(47); %3T
% data{2} = data_extract(51); %4T
% data{3} = data_extract(45); %5T
% data{4} = data_extract(49); %6T

g=figure();
g.Color = Bcolor;
hold on

for i=1:4
    %normalize x with lB which is 10.47 at 6T
    plot_x = data{i}.az_r*((i+2)/6)^0.5/data{i}.lb;
    n_125 = sum(plot_x<125*((i+2)/6)^0.5/data{i}.lb);
    
    %normalized phi with Ec
    plot_y = (data{i}.az_phi-data{i}.avg)*1000/(((i+2)/6)^0.5*data{i}.EC/2.5);
    
    plot(plot_x(1:n_125),plot_y(1:n_125),'LineWidth',0.5,'DisplayName',[num2str(i+2) 'T']);

end
legend();


ax = gca;
box on
ax.Color=Bcolor;
ax.XColor=Fcolor;
ax.YColor=Fcolor;
xl=xlabel('r/l_B');
yl=ylabel('/Ec');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;
ax.YAxis.Exponent = -2;
% axis([0,12.6,-0.6,0.6]);


end
